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.>ESTRji.Cl’ 



This paper reports the results oi an s,nal;>oi.,. oi 
a current repu-latiiiij; cysteii; f03? s. cpcloticon bendj ng 
rr^agnct c.t the 9d’' cyclotron, Lav/rence Radiation 
La'bore.t03:‘5’' , livcrncre, 'California. 

The analj'sis v/as conducted in order to deter:;, i;.e 
v;hy the s;>'stcr' did not regulate sufficiently at certaiii 
times . 

The observations made, find nata obtained i.'.adc it 
.' Ossi'r’le to deterraine the causes f tv'o fac'jors vMch 
did not f.ppo r i. the circuit rep:resentation used in 
design, thoieby reducing the problem to one of correctin 

those t’' o factors namely , noise frora the ■•ener;..tOja , 

-aid tiie operotion of the 1. 'J. amrlif icatio:i at, ^'O, K' 2 * 
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1 . Introduction 



'ihe >G" Cyclotron at the Lawrence Radiation Labora- 
tory, Livermore, California employs tv/o bending magnets 
v/liich are used in certain experiments to direct a beam of 
particles onto a given target. 

Since the target distances are great and the sise 
of the targets small, it is necessa3:y to apply an ex- 
tremely v;ell regulated current supply to the bending 
magnets in order to Jieep the beams on target. 

Lor this purpose, tv;o identical regulating systems 
are used v/ith two motor-generato.r sets to supply the 
current. The generators are rated at 20 KVJ, 160 amps at 
125 volts d.c. The systems v;ere designed to p.i ovide 
stable operation with better than 0 . 1 .% regulated current 
over the range from Ip amps to ICO amps . A diagram of 
the regulating system is Fig. 1-1. An analysis of this 
system v/as to bo carx'ied out, in an attempt to determine 
v/hy it was not performing in a fully a;;ce]>table v;ay . 

The regulator's operation v.’as not . ompletely satis- 
factory and it v/as suggested that the reason might be 
nonlinear operation not accounted for in the original 
design. The original design procedure involved frequency 
response methods v/hich presented no investigation of the 
possibility that nonlinear effects might be important 
over the range of operation. 
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1.0 s'j^steri&uic record oi' doj-icionoi ..c or c.iiu ■ 
ave.ilal';lc , but .ore coran outs rror, c^olot.ron oyoretorb 
u.:!Ci- teclini ciaiic indicated tl'at conietires tire beara 
off target and oscillated., after l^.r'';c charges in ir:..''nct 
current ox^orating level . 

dliG xrlan of attack for tlie i.robloa v/as : 

1. To make a ‘orcliminar'y analj-'sis of the sister.., 
using the original desi.^’rn data, 

2. To ob..rerve the sj-stcja in actual oneration and 
take Gxperinental noasu.rencnts as nccessar^y , 
including i..n/estigation of nonlinoai-itg to 
make certaiii that a suitable cnal;^tical model 
uas established. 

p. To analyze the systera on the basis of that 

model, and ieake recommendations for iraprovcr.mt . 

Since the magnet vras the chief candidate for the 
suspected nonlinearitg , the greatest part of th.e pi-oblem 
becaiae the analysis of the inagnot in order to determine 
the extent to v/hich it v/as nonlinear, ai;d hov; it \.'as 
effecting the o; -ration. The next section exa'lai'xs the 
method of attack emriloyed, and discusses the theory 
behind it . 
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2 . Theory 



In the study of linear systems, the use of the 
block diagram and transfer functions to represent the 
system components is extremely useful, then one of the 
components is a nonlinear olcment, the block diagram 
i/;ay still be used, realisinp; that the transfer function 
concept is no longer valid for the nonlinear block. 

Since such a nonlineai' block v/as suspected to ici.nc 
in the system under study , the attack used in the problem 
included nonlinear techniques v/hich v/ould 'provide a .suit- 
able representation for the nonli'ncari.t;;, , v;hich could be 
used in the block diagram. 

Cne solution is to use the describing function, or 
approximate transfer function v/hich may he defined as 
follov/s. 

Assuming that the input to the nonlinear bloC‘. is 
of the form: 

Input = A sinojt 

then the Courier series fun amental of the output is 

]?(A,(ju) sin ^cut +5 z5(A,o>)J . 

And the describing function is defined as: 




If the approximation is close enough, G^(A,cn) may be 
used in the block diagram. 

A method of evaluating the describing function 
involves taking frequency responses of the nonlinea.i; 
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"block, and using an overlap/ technique with the I'ichols 
chart, to riot a curve of . 

If, in the system under study, the magnet were 
nonlinear, and if it could be represented by a describing 
function, then the system block diagram, could be 

reduced to that shown belov.’, v/here G^^ represents the 
magnet describing function and G-j , the remaining, linear 
part of the system. 




Pig, 2-1 Nonlinear System Block Diagram 

e/? ■ / 

At the limit of stability: 

G"! Gp = -1 



or; 



°1 = Gj 



So that for the amplitude sensitive nonlinearity, 
v/here - is a sing/le curve , stability may be exai. ined 



D 

by checking for intersections of curves of G, and - 

X Wp 

on the Nichols plot. These mig:ht appear as sketched on 
the next page . 
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f/>as<s. 



/ 

Fig. 2-2 Stability Investigation, Using a 
Describing Function 

1 

Fach position along - ^ represents magnitude and 



D 



phase for a particular value of' input marnitude to the 
nonlinear block. 



1. Measure a frequency response for G, , Gj^ at a 
certain amplitude of inj;ut, x, , to'^the non- 
linearity, and make a smooth '^Dode p3ob . 

2. Flot t>"is curve on a !-ichols ch-art, noting 
values of LO (frequency) and x. 

5* -Kepeat for various values of x. 

V/hen G^ 0-^ is plotted on the hichols chart, the 
origin of the hichols plot represents the point 1 /ItO ^ , 
and if G^ Gj^ passes through this point then. 



indicating the liiuit of stability . This point raisht be 
called the critical point. 

If the G^ bloci- is nonlix.ear, the curve G^ G^^ 
v;ill r,,ove v;ith changes in x- 

6 




curve 



Gf G^ = -1 or G^ 



1 



But since G^ is the only 



I 

I 

I 



part v/j'iich changes v.’ith , then f e . ovemei'it of the curve 
must he caused hj 

If instead, of looking at the locus of Gf G^ curves , 
the locus of the critical point relative to the succession 
of G^ Gjj curves is obsei-ved, this locus of critical i;oints 



v;ill be the 




curve . 



kor e::anple : 

ouppose that the sketch bolov; shov/s 
of G^ Gj^ (solid curve) for x = x^, and a 
is at the origin; 



a Kichols plot 



point of 





Pig. 2-p Graphical hetermination at a 
Describing Punction 

The dashed curve vidr-bt represent G^ G^ for x = X 2 * 

llace an overlay sheet O'. er the ’ ichols ■>',iot and trace 

the 0 db and 18C° axes (intersecting at the critical 

point, point C) and the x. curve (including frequency 

7 



locations^. Then adju;..t o' e o , e3:la3' until oils traced 
curve is matched uith the X 2 curve. Kov; mark the origiin 
of the Nichols plot on the overlay ap;air . Tlis new point 
on the overlay sl'ieet represents a second point or the 




?ig. 2-d Graphical Deterrrination of a 
Describing Function 

Gontimie for other values of x until the entire - 
curve is obtained. 

•oo: G. For each succeeding value of x, use the 

overlay technique to plot a corresponding 
point on the _ curve . 

This describing function theory pi-ovided the basis 
for the e;n;erin!cntal './o: .. f'one in section four, which 
v/as airected tov/ai’d evaluating a describing function to 
represent the magnet in the SvStei:: block tia; ram. 
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Chapter I'orr of refei'once tv;o, imalysis ano x^esij?n of 
honlinec-r Feedback Control h> stems Thaler and lactel 
presents a detailed discussion of the describing function 
and its application. 
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3 . 



1 reliridnary .nal;) sis 



The preliminary analysis wo.s carried out usin-y rhe 
results of the original design v/orh, v.'ith some nodii'ica- 
tions as explained below. 



the regulating system along with transf:.:e functions 
determined from the ori'^'inal design v;ork. 

note that the generator and ncejnet each have two 
transfer functions indicated, which are an attempt to 



of the system. Their ori^-^inal design transfer functions 
v!-ere presented as: 



for the magnet . 

These exponents, n, v;ere introduced in an attempt 
to account for the fact that the experimental frequency 
resr.onse curves used to obtain these transfer functions 
fell off to'ward asymptotes with slopes other than 
2C db/decade. 



The block diagiv^m sbov/n in Tic:, q-1 illustrates 



establish tv'o limiting values for the transfer function 




for the generator, and: 




) 
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t « c/c. 



•I'ig. 5-2 gives a cop;>- of the ori'^inal i..a nitiTde 
plot foi- tho genei'ator . Instead of av-corti 'C tl'O tr-n:..for 
functions 




Pig. 5-2 determination of Limiting L'irne Constants 
for the Genex-ator 

\/ith the exponents, tv/o liraiting values for v/ere 

5 5 

obtained v;hich, it v;as hoped, v;ould bracket the actual 
condition. Cne of these >;as gotten Just by setta.ng n=l 
and acceptiiig th.e original value of T^ = (0 .4f) . Then 
^3^5 “ . For the other limiting value, 
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th-G iiiG-Tiitude plot I'orced to approa-.h a 20 db/decacle 

asynptotG a;.- sbo\;n in 5-2, defining the point 

oj 'Tf ' Tf' — 1 

- <^.Z7^ 



from which: 



'T/-- J- 



'■f ~ Zm// ~ZrrCo.S-f) 



4 



SO that: ^ 

Bw the same technique, limiting values of K,,,G,. for the 

*-r 

magnet were obtained as : 

.012 ^ _ .012 

" S (4.0) H- 1 ^^>4^ “ S~l .0795) +'l 



Looking at tPie block 'liagram of ~'ig. 3-1 the re- 
sulting irmiir of transfer functions for the uncompensated 
system became: 

Go = (KiGi)(K 2 )(KjG 3 )(K^G^) 

7.17x10^ 

" Ts+15". 2TCS+2 . 22;TsT.257 

Go’^ = (Kg ) (K 3 G 3 ) * (K.J.G ,J 

5 .88x10^^ 

" T 8 +I 5 . 2 I (S-i-5.64) CS+i2T6j 



The root locus plots for these are' shown in kig. 3-3 
a^vl either T'lot incl-icates that the system is very unstable, 
as expected. 

Introdv.cing the feedback compensation loop, Ak, the 

transfer functions for the compensated system become: 

p 7.18 X 10^' (S+ 9 . 1 )^ 

" (S +13 . ^ ) ( S +499 . 7) (S+ . 08 ) (O' -^4.38718+7257 



and 



p (5.89 X 10^) (S5+9.1)^ 

" TS + 13 . 2 ) ( S+(306 . 7 ) ( S+4 . 56 ) ( S+ ' ■08'5TTS+17T^ 

The regulation desired v;as 0.1%. The s;>stera functions 
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J 



, so thcj.t for the un-starred 



become: 



’ iH-tr 



quantities, ax:>plication of the final value theorem 



leads to the 



fact that 



®7 



(output) 



For 



Jo 



®7 “ 1.CC0C8 



®R 

I.OGOIC 



Both pass the , 1 % regulation requirement by a factor 
of about ten. 

'fhe root loci for. G-^ and are shovm in Fig. 5-4 
and 5- 5 . The G^. locus shov/s the systoia still unstable 
for the gain used. The G^^^ locus iridicatos a stable 
system, but still ve: y lig;htly damped. However, if the 



gain v/ere reduced by a factor of ten, the . 1 % regulation 



specification v/ould still be met, and the system would 
be well enough damped for satisfactory opera.tior . 



Thus if it v/ere acceptable to assume that the actual 
opex'ation v/as bracketed by i the G^, and G-^^ conditions, 
there v/ould arise some doubt as to v/hether the operation 
would bo satisfactory, depending upon v;hetlier G^, or G^,^ 
v/as more nearly correct. 



Bode diagrams v/ere also r/lotted for G^, and G^^^, 
the plot from the original design v/as available for 
comparison . 

In frequency response form: 



since 




) 
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I 



and: 



^ +!)( i ^^ r ^<) (14 ^0 ( i ^, ^0 (1 



/c2^'S- 






These v;ere v;ritten in teicns of i, frequency, instead 
of the usual 6*-> , for coKiparison v/ith the orircinal 

frequency reopjonse plot. 



The magnitude plots for Gj, and compared v;ith 

that presented with the original design, G(^), appear in 
Fig. 3-6. 



The G-r, and G-p* i:-lots do br-acket the G(s) curve , but 
forQ(p), a phase margin of + 50 ° v/as indicated, while those 
predicted by G^, and Gp* were -6.7° and +2° respective!;;,,. 
These results were consi:?tent v;ith the root locus results, 
but conflicted v;ith the design report, casting some doubt 
as to ' liethor the system v/as adequately comp;ensated. 

Flans for the remainder of the analysis v;ere to 
observe aid assess the actus^l oi^eration of the regulating 
system, and to investigate the nonlinearity^ of the magnet 
circuit so that recomm.endations for improvement could 
be made . 
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■T. jXperiiT.ental -jaljsis 

JiOoOd loo; operctiol). ox the v;a. o> .,ervod, 

a) d over i:oat of the operatinf'; ran":e of :: a<^-nct cu.:. rents , 
the system appearod to X^e satisfactory, stable, a.n^ veil 
re-ulated. I.o’«vever, vdo.en adjusted for operation at lov; 
ina net cuiTont, about 20 a.ups , it v;as vossj.hle to ur oduce 
an unsatisfactory condition vrith an osoillatcry ort;;Ut 
and insufficient re ulution ‘by chan;'ir''' the ref ere;-, ce 
level to aive a n.a'"not curorent of al.out ICO .j-yeso 

In the existing sy^sto:: , tv'o filters appear \;hich 
v;ere added fO ‘‘ ad.':> t;i o'cal noise .nav.pr'cssion , a. iter- the 
ori 'inal design v/as coinjjlcted. 

One of these is a r)arallel tee filter \;'i;ich pears 
in the ninor feedhach (compensator) loO'"'? i'- the clock 
labeled in Fig. l-l. ipe filter is locc.tod at the 
input to the jumtode labeled V-h of this amplifier. 

Analysis of this parallel tee, shov/ed that its effect 
is negligible in the operating: range, and can be neglected 
in the system block diagra.’ . 

I'he second filter Lontioneo. is a loa-'-la : filter 
follov;ed by a. ■ arallel tee. This one is located betveeu 
the chopper 6.emo -ulator a, ■! the 12AT7 tube, :'..ch guts 
it in the lorv.'ard path, ri ht ar'ter is; Air. 1-1. 

-1 frequency resTOvse for- this filter vxu recoreed 
aril.- the resiO.tr. indicated t'la.t the phase shift v.’as .. o.j. i- 
gible over the oj)erat:iiig raj ge , but thi-t tie filter .eocuced 
an attenuatioi: of ..hout ton db over the ran -c . 
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x'l'xs attemiation grccxtlj ceci'caa^s tl'C -i i., o_l' 
saj.c uj o / or \'re nt stat ;- acciirao^ requircae t, e.s 

calcul.,.ted in section three. 

In ordex' to invest! 'ate the perfornance ol the 
generator ane espeoiallj- the magnet, by jTrequenc;y lesponse 



techniques , the testing and measuring equii ment \;as set 



up as shovm schematicc^llg in ‘^ig. 1. 

Ihe ..;ost converdent point i'or irtroducin"^’ the test 
si: iic.ls \/as c.t t' c input to F 2 , the r.3. amplifier. 

Jlnce this point happened to be at a lotonti.l or al'out -llC 
volts . uring operatic. the battery and potentiorvT'tor , Idj , 



v/ere introduced to eli:uinc5.te shocji hararc. at the sivnal 
generator, llie .jiynal at t'ue output texninalrj 0 . 1 .' the 
generator v/as extremely noisy - so much so that the noise 
soi.^etiiiies saturated the c-.i'-pli i,'ic:" of the recorder, gl ■ ing 
inaccurate reaeiigs. fhe h-0 filter shown io the diagram 
as inte'o -uced to or.ii!:inate this iiroblein. 

'.Vhe sir^vial e^, from tho ramnet shunt, v/a.: of such a 
S!-:all na.'V.nitude , that the additioial .e.G. amplifier, hj., 
v;as introo.uced to ; laLe e^ readable at the recorder. In 
additioi. , with tho lar,':e amplification reciuirod before 
the rocoxder, it was necessa: y to remove the s all ar.ount 
of d.G. bias at the raagnet shunt, by means of a si. all 
l)attc.ry and potentiometer 
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Frequency' responses for the preneratorr v;ere ta. en 
v.'ith the nia •net disoonnected , and al;. o \:ith' it connected, 
c.s in nornial operation. For both conditions, res^; onscs 
v;ere tahen for valr.es of input, and foi laapnet operating 
levels ovgj the ranr^e in v;hich the reprlator was used. 

I'he generator frequency responses v.'ere ecsentiallj the 
same over the entire operating range, Ihe response for 
the unloaded generator is shovm in 4-2. rig. 4-5 

shov;s that the effect of connecting the marnot is 0 '>^.st 
to reduce the generator volta.ge gain by about six db, 
leaving the form of the resronoo very nearly the same 
as at iio load . Responses for some other conditions a.re 
shown in Appendix I f or ’compa' ison . 

“fitting asymptotes to the loaded generator res/ponse, 
it is seen that the generator may be accur. t^'lp described 
bp the transfer function: 



1 . 1 ‘'-I- . 12 

" (:h75".7??I) = 

Ihe original transfer fu.nction was: 

n ^1- « 8 . 52 ^ 

^5^5 ■ (b/2 .22-fl)‘^ ” (S+2.22)-’ 

For tho mavnet frequency responses, for reasons 
explained in section 2, it \eas desired to take data over 
t])e range of possible magnet currents and. at each of these 
values of 'lagnet euxront, to record the response for 
several different ariiplj tudes of input si.ne v/;_ ve. 

,ith the switch in position one in .-dg. '--1, the sin:. 
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v/ave goneratox' v/aa a^juatcd until the desired, euti.ut at 
the a. gene:cauO,i. terminals (a; o thus the desired ir.Tvt 
at the ..iax'net) v/as observed on the recorder, -.hen c'ith 
the Gi/itci' in ;'o.;ition 2, ma.^'net input and outi'ut rere 
recorded . 



tor e or the many raa^jnet responses ] lotted are 
contained in Appendix II. Froir. these freouency resv onse 
^urn'OS, no consistent varicition could he observed v/hich 
n’ould allov? the determinatioi: o.l a describiiip; function 
by the methods of section tv;o A composite y lot of 

rnapnct rosi-onscs i,_ ii. lig. Bj corapariso;. of the 

results, it \.'as coxicluded that v;ithin e::.perih.ental 
accuracy , the ma iiet fj?enuency res' onse could be con- 
sidered not to h-e am- litude sensitive. Iherefore , an 
a’verage fre ;ucncy resr onso , ■i'". d-p, vas obtained v;hich 

could be used, to X'eprcsont the ns' net for ana-lpsis. By 
fitting acyn'iptotes sr shorn in vh.' . A- C , a transfer 



function may he W itten to approx'- ■ lat the na -net as: 

r _ .0 1C3 (S/2.7- H)(3 /30 -i-1) 

'"m " W.96;i7X’S/T0.7+1T' 

Usinr; these new transfer .fui.iCtion:- the co.m/ ensatedt 
syster.] root locus is Fis . ''Ig corresponding; node 

diapran is '.’i:;';. d— G. Both of these lots ind.ici te that 
the system, should be very stable, and for steady sta'd'^ 
condition, accou,atinr for the attenuation due to 'blie 



filter follov/in'>: B]_^l" 






'R 



T i.coi 



. I e s t e a dy s t a t c re qu i r n: e n t . 
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v.'hich is ri”lt at 
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fro^ the o‘!'_-erv;;t j on o^. the eqv.i; 't i' o "x . , it 

■va^ ioui!’'' ti t :-\:iiioicnt ■■;ain ao.ju. tri. .• u 'c av£ ‘ tic 
at the chopj.er ai;pli±iei to offset the attonai viio; of 
this eatra filtex . 
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'!onclus j.oris 



The e:<:pcrimental analysis conducted or. the generator 
and rnaanet circuit established accurate representations 
for these tu'o i^-rts of the rec^ulatinp: s^/ste’’ . It indicated 
that the maynet circuit v;as not siyniiicantlj an ’■ litu.de 
sensitive in the operating ranyo , and a. linear approxi- 
mation vas obtained. 

Usiny the ne\;' ii.arnet ana generator data, a:^,al5sis 
still indicated that the systen rest onse should be 
satisfactory' , 

-ii.y difficulties in opojration ha to he the result 
of some factor’ not accounted xor in the b'locL cia'ra:;.,. 

'ihe jna.Tict and p'eneraton.’' were cheered experimentally, 
ihe minor feedbach (conpensatin.y) loop resjonse ’..'as also 
-erified experiruentally and aprecd v/cll xith the transfer 
function usee; , ever though an additional ;oa.rallel tee 
filter was in the ci‘:'cuit to reduce noise, ihero 'aa:: no 
reason to doubt the vali-dity of the c’'’opp-''-r amplifier 
transfer function. 

This left tv;o factors which it v;as .felt could have 
caused the trouble-— noise — or the o.G. amplifier, l 2 » 
or both. 

In the desi''"n report, reference one, mention •■■ar; raade 
of a lev; freque.-jo;>^ ripple (about one cps) in the ma'^not 
current, which wa;;. ..ttributed to noise fror t-c ffiTU' ator 
at a f requeue;) of about 59 cps, which idxed with the chopped 
error signal. 'I’he oscillation i^antioned on the first 

pa:'e of section four was of t'‘is form, so tru.t this noise 
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problem na^, not ha’ve been completely elj.minatccl. 

It v;a., ,-ound that the L). ;. am liiicv, 2_-> 
be coroidered as a fixed gain, b.d.-' , over the enti' e 
oporatix.r; rau'-o , as shov/n in the bloc, c.iagra. . I'o 
keep 1,2 constant for all values of ref' rcnce voltage, 
it v;as ne..essar„, to h.al.e a bias adyustnont at ll2» 

Thus, the solution to the problem aj pe-; ’".- to lie in 
moj?e offoctive ;'.'edr:ction .'f noise effects from the •■/•enerator , 
and rei-lucemcnt of Kp by an ampliiri'.r v/hich is effective 
without afyu 'Cl G nij ^ over the entire o"cratiri"' rar^'c of 
reference voltages. 
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Appendix I 



Generaoor Frequency espouse 
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ihe follovrii',; frequency reS; onses aie for the 
generator at difxerent voltage outpute . 'I'ho response 
v;as noarl.y tlie same for all output voltages. 

All responses v.’ore taken with inputs at the d.c. 
amplifier, F 2 • The 6C volt curve, page 39, includes 
a 3-7 ‘■'-’b gain due to I' 2 * The 3^' volt curve, pi.gc 38, 
has been coriocted, hy subtracting K 2 (dh) . 
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Appendix II 

-lagnet /'requency Kenponses 



^■0 



The follo'./ing curves are some frequenc;)' resro.ises 
for the raagnet, at cifferent operating conditions did 
input amplitudes. These are so;ne of the responses used 
in the rnapn' t analysis of section four*. Cn each curve 
is indicated the mavnet current level for the response, 
and the value of constant input in volts, maintained 
at the input to the magnet throughout the run. 
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